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t'r'e  w,  r'r  throegi;  30  September  1974  was  submitted  In  Deccinher  of  1974.  Con- 
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N'lrOl  -0216-0033  and  during  the  period  of  .N00014-73-C-1]  78. 
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The  initial  phase  of  an  experimental  data-taking  program  has  been  con- 
pleted  lending  to  the  measurement  of  the  spatial  correlation  function  of  the 
phase  of  an  .'IF  wave  after  passage  through  the  ionosphere.  Tlie  propagation  path 
for  tr.e  experiment  was  West-East,  a distance  of  2,480  kilometers.  The  transmit- 
ting site  was  in  Boulder,  Colorado.  The  receiving  site,  which  was  in  the  Dela- 
ware Valley,  Pennsylvania,  had  a 40-kilometer  baseline.  The  receiving  system 

consisted  of  8 phase-synchronized  HF  receivers.  Hie  raw  data,  which  are  in  ex- 
g 

cess  o:  10  bits,  have  been  recorded  digitally  on  magnetic  tape.  Data  reduction 
and  analysis  will  take  about  one  year. 


IirTRODL'CTION 

Nonhonogeneities  in  the  refractive  index  of  the  ionosphere  induce  time- 
varying  spatial  variations  in  the  phase  of  a wave  received  from  a distant  high 
frequency  (HF)  radio  source  after  passage  of  the  wave  through  the  ionosphere. 
TT.e  correlation  distance  of  the  phase  approximately  equals  the  nominal,  maximum 
useful  aperture  length  of  an  HF  antenna  whether  it  is  designed  for  point-to- 
point  communications  or  for  over-the-horizon  radar.  In  the  former  case  large 
antenna  size  helps  in  providing  a large  s 1 gnal-to-nolse  ratio  (S"K)  which  In 
turn  is  required  for  the  maintenance  of  a useful  communications  circuit.  When 
the  -rir.'.t :.:ia  size  is  less  than  the  nominal  correlation  distance,  the  SHR  grow.s 
linearly  '..'ith  antenna  size;  this  is  because  all  signal  components  reccivod  by 
t'l--  antenna  are  surcim-d  coherently.  Howi‘ver  when  the  antenna  size  exceeds  the 
nominal  correlation  distance,  signal  decorrelation  across  tl»e  aperture  pre- 
cludes fullv  coherent  sunuiuition  and  re.sults  in  less  than  a linear  buildup  of 
SNR  wit;,  antenna  size. 


A 
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In  the  sjocuntl  case,  that  of  over-the-hor  Lzon  radar,  a .secon<I  factor  of 
import. ’.nee  arise;..  In  addition  to  a linear  improvement  in  SI;r  with  antenna 
si.'e.  the  angular  directivity  similarly  lt:iproves  (the  beamwidth  of  a diffrac- 
tion limited  aperture  is  approximately  the  ratio  of  wavelength  to  aperture 
site).  However,  wh.er.  t'ne  antenna  length  exceeds  the  spatial  correlation  dis- 
tance, ti'.e  aperture  ceases  to  be  diffraction  limited  and  the  angular  direc- 
tivity no  Icnger  grov.-s  with  size.  It  is  in  tiic  interest  of  overcoming  some 
rt  t'.e  limitations  in  ti.e  i!F  radar  problem  that  an  experimental  40-kilometer 
iH'"  cata-taking  array  v.as  constructed 

'.’here  is  a third  case  in  which  knowledge  of  the  spatial  correlation  func- 
tion may  prove  useful  in  future  system  design:  It  is  in  the  field  of  HF  direc- 
tion finding  (’HFDF),  Tlie  apparent  bearing  of  an  HF  source,  as  determined  by 
the  direction  of  arrival  of  an  liF  wave  from  that  source,  often  appears  to  wan- 
der within  a few  degrees  of  the  true  bearing.  Hie  source  of  the  apparent  bear- 
ing wander  is  large-scale  wave  motion  within  the  ionosphere,  sometimes  called 
travelling  Ionospheric  disturbances.  'ihe  long  wavelengths  (tens  of  kilometers) 
of  these  internal  waves  in  the  ionosphere  and  their  speeds  (the  order  of  100 
meters  per  second)  result  in  a time  scale  of  bearing  wander  which  is  the  order 
of  minutes.  Observation  of  the  apparent  bearing  over  several  such  cycles  tends 
to  disclose  the  true  bearing,  in  essence  by  averaging  the  bearing  data  over 
t im.e . 

;'u»  average  of  bearing  data  also  may  be  r.ade  over  space,  provided  that  the 
r.casurem.ents  extend  over  a suitably  large  distance;  that  is,  over  a measuring 
interval  large  compared  to  the  spatial  correlation  distance.  Hence  a nondif- 
fr act  ion- limited  aperture  may  prove  useful  in  HFDF.  /Mthough  an  optimum  signal 
processing  theory  for  large  aperture  HFDr'  has  yet  to  be  developed,  one  of  tlie 
dominant  par.ameters  of  the  theory  is  known  to  be  the  spatial  correlation  dis- 
tance of  the  phase  of  the  HF  wave  after  it  pas.ses  tlirough  the  ionosphere.  Thus 
the  information  to  be  extracted  from  the  data  may  prove  useful  in  HFDF. 

The  Valley  Forge  Research  Center  of  the  Moore  School  of  Electrical  Engin- 
eering, Iniversity  of  Pennsylvania,  is  engaged  in  an  ongoing  study  of  the  lim- 
itin’, angular  resolution  properties  of  large  apertures,  and  of  self-correction 
techniques  for  overcoming  the  effects  of  nondif fraction-limited  and  time-varying 
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nperturfls  fl].  In  ii  two-year  stndy  of  the  IIF  radar  problem,  a theory  was  de- 
veloped of  the  naxiiTum  aperture  wliich  could  be  made  to  function  as  If  it  were 
diffraction  limited  throup.h  the  use  of  adaptive,  self-focusing  signal  proces- 
sing techniques  subject  to  the  conditions  experienced  by  IIF  signals  in  passage 
through  th.e  ionosphere  [2,3].  Tlie  Raximum  size  was  found  to  be  proportional 
to  the  phase  correlation  distance.  Ivlien  a search  of  the  literature  disclosed 
that  t'.'.is  quantity  had  hitherto  not  been  measured  (perhaps  because  of  the  in- 
herent difficulties  in  the  measurement)  an  experimental  system  was  developed 
to  measure  the  necessary  data,  and  five  successful  data  recording  experiments 
have  been  conducted.  Tliese  data  are  stored  digitally  on  magnetic  tape.  Later 
we  expect  to  derive  from  these  data  the  relevant  spatial  statistics  of  the 
phase  of  the  HF  signal  received  via  the  ionosphere  during  these  five  experi- 
ments. 

THE  EXPERIMENT* 

An  ionosonde  of  the  Environmental  Research  Laboratory,  National  Oceanic 
and  Atmospheric  Administration,  U.  S.  Department  of  Commerce,  Boulder,  Colo- 
rado, radiated  an  HF  pulse  train.  The  carrier  freq".ency  was  9.010  MHz,  the 
pulse  repetition  frequency  was  lOOHz , and  the  duration  was  64  microseconds. 

Both  the  carrier  frequency  and  the  repetition  frequency  were  deriv'ed  from  an 
atomic  standard  and  hence  were  highly  stable.  The  transmitting  antenna  was 
modified  so  as  to  depress  the  beam  from  the  vertical  in  an  easterly  direction, 
thereby  permitting  a signal  with  adequate  signal  strength  to  be  received  in 
the  Delaware  Valley,  Pennsylvania. 

Hie  receiving  array  consisted  of  eight  receivers  in  the  Delaware  Valley. 

(1]  See  articles  1 and  2 of  V^FRC  (Quarterly  Progress  Reports  Nos.  8 through  14. 

(2]  S.  Hassan  Taherl,  "imaging  through  Turbulent  Distorting  Media",  VFRC  QPR 
Nci.  9,  pp.  47-65,  .June  1974. 

(3]  S.  Hassan  Taherl,  "Experiment  for  the  Measurement  of  the  Spatial  Struc- 
ture Function  of  the  Phase  Fluctuations  due  to  Ionospheric  Propagation", 

VFRf;  OPR  No.  10,  pp.  36-47,  August  1974. 

* 

Details  on  the  experimental  equipment  arc  found  in  Valley  Forge  Research 

Cienler  Quarterly  Progress  Reports  Nos.  10  through  14,  (August  1974  - August  1975). 


Tlu*  recoiviTS  wen’  special  units  desiyiied  and  built  by  VI'KC  staff  for  this  ex- 
peri~e;il  . Tlie  naxirtun  separation  was  nearly  AO  kilometers.  The  locations  and 
spacirots  ’■•■ere  so  chosen  that  th.e  set  of  inter-pair  spacings  would  provide  a use- 
ful distribution  of  abscissa  values  for  sample  points  on  tiie  phase  correlation 
function.  Ike  positions  are  given  in  the  map  of  Figure  1.  Table  1 lists  the 
distances  (km)  and  Table  2 lists  the  directions  (degrees  clockwise  from  North) 
between  stations.  The  locations  were  at  the  Research  Laboratory,  at  homes  of 
: ic’ulry  and  staff  members,  at  the  University's  Flo’wer  and  Cook  Obser’/atory , and 
at  a nearby  brass  foundry. 
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FIGURE  1.  HF  RECEIVING  ARRAY 
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i';;e  funcLion  of  the  rcceivinp,  system  was  to  measure  and  to  record  phases 
of  I’.F  signals  from  Boulder  received  by  the  eipht  receivers.  Accomplishment 
o:  th.is  task  required  synchronizing  the  local  oscillators  in  the  receivers, 
r-.ese  rh.ase  synchronized  reference  waves  were  obtained  from  the  carrier  of  a 
local  radio  station,  KYV,  which  broadcasts  an  acceptable  signal  strength  at 
1160  :j1z  throughout  the  Delaware  Valley.  Within  each  receiver  the  carrier 
was  freed  of  its  modulation  by  a crystal  filter,  resulting  in  a 1.06  MHz 
reference  wave  with  appro.ximuitely  1°  rms  phase  jitter.  A frequency  synthe- 
sizer in  each  receiver  developed  a large  set  of  stable  frequencies  in  the  HF 
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region  of  the  spectrum  from  this  phase  stable  reference  wave.  Proper  choice 
of  transmitter  frequency  and  local  oscillator  receiver  frequency  permitted 
phase  detection  of  the  received  wave. 

Figure  2 shows  the  organ ijiat ion  of  the  receiving  system.  A phase 
ccnparisan  is  made  between  the  local  reference  and  the  succession  of  HF 
pulses  iron.  Colorado.  Hie  analog  phase  measurement  in  each  receiver  was 
sampled  iOO  times  per  second,  digitized,  and  transmitted  via  data  modem 
and  telephone  circuit  to  the  Valley  Forge  Laboratory,  where  the  signals  were 
entered  into  a buffer  digital  store,  adjusted  in  format,  and  read  out  onto 
a 7- track  tape. 


FIGURE  3.  RECEIVING  INSTALLATION  AT  SITE  8 


Figure  3 shows 
scope  photo  showing 
timing  marker  (just 
tracking  gate. 


the  receiving  installation  at  Site  8.  Figure  4 is  a 
the  envelope  of  tlie  received  pulse  (upper  trace)  and  a 
under  the  main  pulse)  showing  tlie  location  of  the  pulse 
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FIGURE  4.  OSCILLOSCOPE  PHOTO  SHOWING  RECEIVED 

SIGNAL  (UPPER  TRACE)  AND  TRACKING  GATE 


Five  successful  tests  were  run  during  June,  1975.  Tne  dates,  times 
and  durations  are  given  in  Table  3.  (Midpath  time  is  approximately  one  hour 
earlier,  and  GMT  is  five  hours  later.)  Hie  recorded  data  are  in  excess  of 
10^  bits. 


DATE 
12  June 
21  June 

24  June 

25  June 
27  June 

TABLE  3 


TIIS:  (EST) 
6 : 30  A.'-; 
5 : 30  AM 
4 : 00  AM 
Midnight 


DURATION  (min) 
60 
120 
60 
60 
60 


7 : 00  AM 
SUCCESSFUL  TESTS 


Tne  primary  statistic  sought  in  this  work  is  the  spatial  correlation 
function  of  the  phase  of  the  received  wave  after  its  passage  through  the  iono- 
sphere. The  abscissa  of  the  correlation  function  is  a distance  on  the  ground. 
Samples  of  the  correlation  function  will  be  measured  at  distances  equal  to  the 
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spacings  between  receiving  stations.  Correlation  information  in  the  East-West 
direction  (along  the  direction  of  propagation)  and  in  the  Korth-Sonth  direction 
(transverse  to  the  direction  of  propag.ation)  will  be  obtained.  On  theoretical 
grounds  the  correlation  functions  are  expected  to  be  similar.  The  resources  of 
the  experimental  system  were  devoted  roughly  equally  to  measurements  in  these 
f-o  directions  and  with  som.e  overlap  to  permit  comparison.  Table  4,  construc- 
ted from,  the  earlier  tabular  data,  shows  the  interstation  distances  in  the  K-S 
and  K-W  directions.  Fourteen  other  interstation  separations  permit  measure- 
ment of  correlation  coefficients  in  a variety  of  directions;  comparison  v;ill 
be  made  between  these  values  and  those  measured  along  the  direction  of  propa- 
gation (E-W)  and  transverse  to  it  (N-S). 


N - 

S 

E - W 

STATIONS 

SEPARATIONS (km) 

STATIONS 

SEPARATIONS (km) 

1-2 

2.0 

2-7 

9.9 

3-4 

4.7 

2-6 

11.6 

2-3 

3.9 

7-8 

15.3 

1-3 

5.8 

6-7 

21.5 

2-4 

7.4 

2-8 

25.3 

1-4 

9.3 

6-8 

36.8 

5-7 

27.3 

5-8 

28.1 

TABLE  4.  VALUES  OF  THE  ABSCISSA  (IN  KM)  OF  THE  COORDINATES 
OF  THE  SAMPLES  OF  THE  SPATIAL  CORRELATION  FUNCTION 
OF  THE  PHASE. 


Tlie  value  of  the  correlation  function  at  a distance  corresponding  to 
the  spacing  between  the  ith  and  j th  stations  may  be  calculated  either  by  direct 
multiplication  or  from  the  variance  of  the  phase  difference  function.  The  cor- 
relation coefficient 
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where  arc  the  time  samples  of  the  received  phase  after  removal  of  the 

artifacts  and  the  slowly  varying  pliaso  roll  compy’ nc.nts  ascribable  to  slight 
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fri‘c,uency  dif  fcroni'i-s  In’twoiMi  the  transpitter  and  receivers.  llie  overbar  re- 
presents a tine  averap.e.  Since  the  pliise  variations  and  are  derived  from 
a conmon  source,  it  is  assumed  (in  the  expression  above)  that  their  mean  square 
values  are  the  same;  this  assumption  will  be  tested  by  independent  measurement 


e: 


ror.  (1)  that  p^.  is 
also  is  available 


available  from  the  direct  measurement 
from  the  variance  of  the  phase  differ- 


cr.cc;  i 


'ij  " (2) 

Although  (1)  is  more  direct,  there  may  be  instances  in  which  (2)  will  prove 
useful,  as  when  rapid  movement  in  the  height  or  the  tilt  of  the  ionosphere 
causes  rapid  but  correlated  variations  in  the  phases  of  the  received  waves  at 
two  receivers.  Such  correlated  variations  are  eliminated  in  calculations  based 
on  the  phase  differences. 

Hie  noise  properties  of  the  two  measurement  procedures  are  different,  a 
factor  which  has  to  be  reckoned  with  in  the  numerical  calculations.  If  and 
n^  are  the  time  samples  of  adaptive  phase  noise  in  the  ith  and  j th  receivers, 
the  average  product  of  the  phase  measurements  becomes 


Y 


(y*!  + n^)  + n^) 


= d . i> . + i.n.  + b.n.  + n.n.  = 

1 J i J 1 1 1 J i 3 


(3) 


which  is  unbiased.  The  latter  three  terms  reduce  to  zero  under  the  assumptions 
(1)  that  a received  wave  and  the  receiver  noise  are  Independent,  (2)  that  both 
are  zero  mean  and  (3)  that  the  phase  noises  from  receiver  to  receiver  are  sta- 
tistically independent. 

'Hie  phase  difference  variance,  on  the  other  hand,  is  not  unbiased  by  the 
noise.  Equation  (2)  becomes 

Ojj  = (d^  + - dj  - rij)2  = 2<j>^  - 2dj^dj  + (4) 


under  the  same  assum.pt  Ions . Under  a fourth  assumption  that  the  additive  phase 
noises  In  tlie  various  receivers  share  common  statistics,  (4)  becomes 
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O?.  = 2P  f 2.r  - 2«^«. 


(5) 


The  phase  noise  power  can  be  measured  usinp,  (5)  provided  that  two  of 
the  stations  are  sufficiently  close  to  each  other  such  that  their  phases  have 
nearly  unity  correlation,  in  which  case  ^ also  can  be  measured 

when  two  stations  are  sufficiently  far  apart  that  their  phases  have  nearly 


aero  correlation,  in  which  case  v^'i’ . “ ~ ® under  assumption  (2)  above; 

ccr.paripcn  cf  (3)  and  (5)  then  yields  n'^.  If  both  situations  can  be  found 

to  exist  in  the  data,  then  independent  measurements  of  can  bo  made  and  com- 
pared . 


.•\nother  interesting  statistic  which  may  be  available  from  the  phase  data 
is  the  first  probability  density  function  (pdf)  of  the  signal  phase.  The 
measured  pliaso  in  an  arbitrary  receiver  is  ^ + n where  is  the  signal 

phase  and  n is  the  additive  phase  noise.  Eased  on  the  first  assumption  made 
earlier  (statistical  independence  of  (J  and  n) 


w(5)  = w(';)  * w(n)  (6) 

where  '.K")  means  pdf  and  * means  convolution.  Tfius  the  shape  of  w($),  whicl> 
can  be  extracted  from  the  data,  differs  from  w((l)  , the  pdf  of  the  signal  phase, 
because  the  latter  is  convolved  with  the  pdf  of  the  phase  noise. 

Although  shape  information  nay  be  obscured,  the  signal  phase  variance  can 
be  calculated  from  the  variance  relation  in  summing  independent,  random 
variables  : 

0?  = (7) 

iji  p n 


Based  on  assumption  (2)  above  (zero  ineans  of  p and  n),  a^=  ^ and  = n^, 

T n 

the  calculations  of  which  were  discussed  earlier.  Thus  the  signal  phase  vari- 
ance is  the  difference  between  the  mean  square  value  of  the  total  phase  mea- 
surement and  the  value  calculated  for  the  noise. 

Shape  information  also  may  be  possible  to  extract.  If  n <<  p , as  is 
expected,  w(p)  will  be  slightly  broader  than  w(p)  (according  to  (2)).  Tiieir 
shapes,  however,  will  be  very  similar.  In  particular,  if  the  phase  noise  power 
Is  sriali  compared  to  tlie  signai-pliase  power  and  if  w(p)  appears  Gaussian,  It 
may  be  assumed  that  w(p)  also  is  approximately  Gaussian. 


■f" 
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BEST'AVAIUBIE  CO?Y 

uiscrssroN 

The  spatial  variation  of  t!;e  phase  of  a signal  received  along  a line  on 
tlu'  ground  after  ih.e  wave  has  passed  through  the  ionosphere  is  believed  to  have 
at  least  two  components.  The  first  is  due  to  scattering  from  random  nonhonogen- 
eities  cue  to  rar.dcm  ionization  and  rccocbinatlon  phenomena.  Tlie  second  is  due 
to  large-scale  effects  resulting  from  wave  motion  in  the  ionosphere.  The  former 
introduces  a rancor,  phase  variation  across  an  aperture  which  sets  an  irreducible 
limit  to  the  effective  maximum  diffraction-limited  length  which  may  be  achieved. 
Ti'.e  l.attor  is  quasi-dcterministic  in  a quasi-static  sense;  that  is,  tine-varying 
curvatures  and  tilts  in  the  ionosphere  due  to  the  passage  of  waves  may  induce 
relatively  simple  phase  variations  in  the  wave  front.  For  example,  at  a given 
Instant  in  time  the  phase  variation  of  a wave  across  a large  aperture  night  be 
well  represented  by  a low  order  polynomial.  The  constant  term  represents  the 
fixed  phase.  The  linear  term  represents  an  angular  tilt  of  the  ionosphere.  The 
quadratic  tern  represents  a curvature  associated  with  the  mean  refractive  index 
of  a ray  path  through  the  ionosphere  as  a function  of  transverse  position  with- 
’ in  it.  This  term  nay  be  identified  with  the  properties  of  a spherical  lens. 

High  order  terms  would  be  the  equivalent  of  distortion-inducing  features  of  an 
optical  system. 

The  spatial  correlation  junction  which  we  will  attempt  to  extract  from 
the  recorded  data  will  include  both  the  effects  of  scattering  and  the  large- 
scale  phenomena.  The  result  of  the  measurements  will  set  bounds  or  lower  lim.its 
to  the  maximum  size  that  diffraction-limited  apertures  can  be  made  under  con- 
ditions in  which  the  measurements  were  made.  However,  through  the  use  of  rela- 
tively sophisticated,  adaptive  signal  processing  circuits  of  the  type  being  ex- 
plored at  V'FRC  some  of  the  normally  deleterious  effects  of  the  large  scale  phe- 
nomena may  be  removed.  Specifically  we  believe  that  the  linear  term  may  be  re- 
moved on  a continuous,  dynamic  basis.  Ke  feel  too  that  the  slowly  changing 
quadratic  variation  can  be  removed.  Thus  the  irreducible  distortion  which  would 
limit  useful  size  of  an  aperttirc  would  be  determined  by  the  contributions  to  the 
spatial  correlation  of  the  phase  of  the  energy  scattered  by  the  small-scale, 
random  variations  In  refractive  index  within  the  ionosphere  and  by  the  higher 
order  terms  of  a polynomial  expansion  of  the  ph.ase  variation  induced  by  the 
large  scale  ionospheric  phenomena. 


-li- 


lt is  not  knouTi  at  this  writing  the  extent  to  wiiicii  tlicsc  plienomena  will 
bo  separable  in  the  data.  The  calculation  of  the  spatial  correlation  function, 
as  described  earlier,  ignores  the  distinction  made  in  this  section.  During  the 
course  of  the  data  reduction  program,  consideration  will  be  given  to  the  separa- 
tion of  the  effects  of  these  different  phenomena.  Attention  will  also  be  given 
to  the  development  of  any  unexpected  statistical  or  deterministic  phenomena  ev- 
ident from  the  data. 
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JO  abstract  (Con„nu,  Id.  If  n.c...^  t,r  "“■"'’•'^'^Experimental  equipment  was 

designed,  built,  and  used  in  an  experiment  for  the  purposes  of  acquiring  data, 
performing  analyses,  testing  algorithms,  and  evaluating  system  characteristics 
for  measurement  of  the  spatial  correlation  function  of  the  phase  of  an  HF  wave 
after  passage  through  the  ionosphere.  For  the  test,  8 phase-synchronized  rcceiv 
ers  were  designed  and  built,  distributed  randomly  over  a 40-kilometer  baseline, 
and  v.sed  to  receive  HF  signals  from  a transmitting  site  some  2500  kilometers 
distant.  Raw  data  in  excess  of  10  blt.s  were  . collected^  
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